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ABSTRACT 

We present a new semi-analytic model for the common growth of black holes (BHs) and 
galaxies within a hierarchical Universe. The model is tuned to match the mass function of 
BHs at z = and the luminosity functions of active galactic nuclei (AGNs) at z < 4. We use a 
new observational constraint, which relates the luminosity of AGNs to the star-formation rate 
(SFR) of their host galaxies. We show that this new constraint is important in various aspects: 
a) it indicates that BH accretion events are episodic; b) it favours a scenario in which BH 
accretion is triggered by merger events of all mass ratios; c) it constrains the duration of both 
merger-induced star-bursts and BH accretion events. The model reproduces the observations 
once we assume that only 4 per cent of the merger events trigger BH accretion; BHs accretion 
is not related to secular evolution; and only a few per cent of the mass made in bursts goes into 
the BH. We find that AGNs with low or intermediate luminosity are mostly being triggered 
by minor merger events, in broad agreement with observations. Our model matches various 
observed properties of galaxies, such as the stellar mass function at z < 4 and the clustering of 
galaxies at redshift zero. This allows us to use galaxies as a reliable backbone for BH growth, 
with reasonable estimates for the frequency of merger events. Other modes of BH accretion, 
such as disk-instability events, were not considered here, and should be further examined in 
the future. 
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1 INTRODUCTION 

Most massive galaxies in the loc al Universe host supermassiv e 
black holes (BHs) in their centres jKormendv & Richstone|[l995h . 
Active periods of accretion onto such objects, seen as active galac- 
tic nuclei (AGNs), are observed in less than a few percent of the 
objects. Consequently, models of BH evolution rely on assump- 
tions about the triggering of accretion events and the connection, if 
any, with the host galaxy evolution. A key problem in such studies 
relates to the big difference in scale between the cold gas in the 
galaxy (distributed over several kpc) and the inner disk around the 
BH (typically smaller than a parsec). Various suggestions have been 
made to bridge the gap between these scales, including: merger 
triggered accretion, secular processes in the host galaxy including 
disk instabilities mass loss from evolved stars, star formation (SF) 
activity in the central region, and more. 

The merger scenario is motivated by two different lines of ev- 
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^ In this work we use the terms "secular evolution" and "secular processes" 
to indicate all types of internal processes within a galaxy. These include 
events of disk instability that might induce bursts of star-formation. Secular 
evolution does not include events related or triggered by galaxy mergers. 



idence. First, A^-body and hydrodynamical simulations of merg- 
ing galaxies suggest that such events drive a large amount of cold 
gas into the centre of the merging system, re sulting in the forma- 
tion o f a large bulge or a spheroidal galaxy (Barnes &Hemauist 



F 



tion or a large bulge or a spheroigal galaxy (Itiames & i-Lemguisi 
1991'; Mihos & Hernquist 'l996'; ICox et al.1120061 : [Robertson et al 



2,006; Di Matteo et al. 2007). Second, AGN activity is observed in 
local ultra-luminous infrared galaxies (ULIRG s), which are known 
to be as sociated with ener getic merger events ("S anders et al.lll98 



Sanders & Mirabel 1996; Surace et al . 2000; Ca nalizo & Stockto 



200ll : IVeilleuxetal.ll2009l) . Observationally, an evidence for a di- 



rect connection between minor mergers and AGN activity is more 
difficult to establish. 

Recent studies show that most galaxies at z < 2 that 
host a low-luminosity AGN do not show morpho logical evi- 
dence for major merger events in their recent historv (iGabor et al 



dence tor major merger events m tneir recent mstory (luabor et al 
20091 : ICistemas et alj201ll : iKocevski et aLll2012l : ISchawinski ~ 



20121). Only high-lu minosity AGNs are related to major mergers 
dXreister et al.l2012l) . This has been taken to indicate that the role of 
secular processes in driving gas to the centre of galaxies and feed- 
ing the BH is iniportan ( : and very co mmon at almost all redshifts 
(lEfstathiou et al.' 1982; Genzel et al. 2006: 'El megreen & BurkertI 
[2010 ; Ceveri no et al, 2010; Bourn aud et al. 2012j). The nature and 
frequency of such internal events are not fully understood, partly 
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because it is difficult to discriminate them from minor merg- 
ers. Nonetheless, various structural properties of galaxies have 
been suggested to dist inguish between secular and merger events 
jKormendv et alj2009l) . 

Most of the studies mentioned above have measured the mass 
ratio in mergers by visual inspection, involving a large uncertainty. 
Moreover, there may well be a time-delay between the merger event 
and the peak luminosity of the AGN. Since recently merged galax- 
ies ha ve similar clusterin g properties as other galaxies of the same 
mass teonoli et alj201(3h . the effect of time-delay is difficult to de- 
tect (h owever, it might be constrained by detailed analysis of g alaxy 



spectra lDavies et al.ll2007l : Iwild et al.ll2007l : IWild et al.ll201[ 

Other internal mec hanisms within a galaxy were proposed as a 
trigger for BH accretion. iHopkins & Hemquistl ( l2006l) have demon- 
strated th at feedback from the BH can regulate low levels of BH 
accretion. iBest et aL have argued that observations of radio- 

loud AGNs provide indications for hot gas accretion into BHs. 

The mass of super-massive BHs in galactic centres is st rongly 
correlated with both the bulge mass ("Magorrian et al.' 1998; 



McLure & Dunlotj I2OO2I: iMarconi &liunt 2003 : Hiiring & Rix 
2004ISani et al.ll201 ll: iMc Connell & MH2012h and the bulge ve - 



locity dispersion jperrarese & MerrittlbOOOl: iGebhardt et alj|2000t) 



These relatively tight correlations suggest that the mechanisms 
which are responsible for building up the mass in bulges are also 
related to BH growth and hence the triggering of AGN activity. 
However, since the mass is an integrated quantity, these corre- 
lation s could also or iginate fr om the hierarchical agg regation of 
mass jnirschmann et al. 2010: Jatinke & Maccid[201ll) . Better un- 
derstanding of the various accretion scenarios is therefore required 
in order to understand the BH-bulge relationships. 

Further hints for the nature of BH accretion could be pro- 
vided by relating AGN luminosities to SF rates (SFRs) of their 
host galaxies. Unlike morphological studies, SFRs are well de- 
fined, and can be quantitatively compared between various pop- 
ulations of galaxies. Studies of SFRs and AGNs luminosities in 
large s pectroscopic surve ys, like the Sloan Digital Sky Survey 



within star-formine galaxies (iKauffmann et al. 20031; Wild et al.l 


20071 Salim et alj 200l Daddi et alj |2007|; Da 


ivies etalj 20071; 


Silverman et all 20091; Netzer 20091; IWild et al.l 


2OIOI; iMor et al.l 


2OI2I; Rovilos et al. 20121; Tommasin et alj 


2012). Here we choose 


to highlight the recent measurements by 


Shao et al.l (I2OIOI) and 



iRosario et all l [2bl2l) that are based on SFRs obtained by Herschel 
in the FIR, and are more r eliable than UV-based SFRs. Accord- 
ing to iRosario et al.l | |2012|) . such SFRs are uncorrelated with the 
AGN luminosity at low and intermediate luminosities and at high 
redshift. However, high luminosity AGNs at low and intermediate 
redshifts do show a coiTelation with the SFRs of their host galaxies. 

In this paper we use a semi-analytic model (SAM) of galaxy- 
formation combined with BH evolution, in order to examine the 
congelations between SFRs and AGN luminosities. SAMs are nat- 
urally being used for this problem as they provide a statistical 



I Kauffmann & HaehnelJ 200Cl; Croton et alj 20061; Bower etal. 


200e 




Malbon et alj 20071: iMonaco et al. 


2007: Somerville et al. 


200s 


Bonoh et al. |2009|; iFanidakis et al 


I2OIII; IFanidakis et al. 


2012 


:lHirschmann et all 20121). 



Our study is different from p revious SAMs in various as- 
pec ts. First, we use th e SAM from iNeistein & WeinmannI ( I2OI0I) 
and lWang et alj ( 1201 2|) that fits the stellar mass function of galaxies 
at < z < 4 and the clustering of galaxies at low redshift. As a re- 
sult, the star-formation histories and merger-rates of galaxies within 



the model are in broad agreement with the observed Universe. Sec- 
ond, we do not attempt to model AGN feedback as a function of 
each specific AGN but rather assume that the galaxies within our 
model experience an average AGN feedback that depends only on 
halo mass. Third, we aim to match a wide range of AGN obser- 
vations, including the relation between SFR and AGN luminosity 
mentioned above, and the fraction of host galaxies that are experi- 
encing major merger events. 

Various recent empirical models of BHs and AGNs evo- 
lution are aimed to fit th e BH mass function and AGN 
lumin o sity functions (e.g . Wvithe & Loebl 20031; Lapi et all 
20061; iHopkins et all l200d I200I I2OO8I ; IShankar et alj |2009| 



ShankaJ 1201 Ol IShankar et all 1201 2t IConrov & Whitel I2OI2 : 
Draper & Ballantvnell2012h . Here we choose to use SAMs, as they 
allow a more direct modelling of SFRs and merger fractions of in- 
dividual objects, in close contact with the formation history of their 
host haloes. These features play a key role in our model. 

This paper is organized as follows. In section |2] we describe 
the model ingredients and the details of computing all the proper- 
ties of galaxies and BHs. Section [3] includes a detailed compari- 
son between the model and observations. The results are summa- 
rized and discussed in section |4] This study is based on the cos- 
mological parameters that are used by the Millennium simulation: 
{Q.m,O.A, 0-8, h) = (0.25, 0.75, 0.9, 0.73). Throughout the paper 
we use log to denote logjQ. 



2 THE MODEL 

2.1 Galaxy formation and evolution 



The SAM used in this work is ad opted from | Wang et al.l ( 2012 ). 
and is based on the formalism of INeistein & WeinmannI ( 20IC ). 



The model follows galaxies inside the complex structure of subhalo 
merger-trees taken from a l arge A'^-bo dy cosmological simulation 
(the Millennium simulation. ISpringel e t al. 2005). This simulation 
follows 2160'^ dark-matter particles inside a periodic box of length 
500 /i^^Mpc, with a minimum halo mass of 1.72 X 10^*^ H-^Mq. 

Our SAM includes the effects of cooling, star formation (SF), 
accretion, merging, and feedback. Unlike other SAMs, these laws 
are simplified to be funct ions of only the host su bhalo mass and red- 
shift. We have shown in INeistein et alj ( |20I2|) that this concept is 
sufficient for the SAM to reproduce the gas and stellar mass content 
of galaxies, on an object by object basis, as obtained from a hydro- 
dynamical simulation to an accuracy of 0. 1 dex. Consequently, this 
SAM is complex enough to accurately follow the SF histories of 

galaxies. 

The specific model used here was pres ented in 'Wang et akl 

2OI2I) s model 4. It is based on the SAM of iDe Luci a & Blaizgn 
20071) . transformed to the language of INeistein &^Weinmannl 
20101) and including the following modifications: 



• SF efficiency (the ratio between the SFR and the cold gas mass 
within the disk) is lower in small mass haloes. 

• Satellite galaxies experience stripping of hot gas in propor- 
tion to the stripping o f dark-m atter, following the suggest ion of 
IWeinmann et al T( l2010l) (see also lKhochfar & Ostrikeill2008l) . 

• Cooling is suppressed within haloes that are more massive 
than 5.6 x 10^^ /i~^M0, and there is no SF within haloes of mass 
> 5 X 10^^ /i~^M0 at z < 1.3. These ingredients are aimed to 
mimic the feedback from AGNs, although they do not depend on 
the specific AGN hosted by each galaxy. 
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• The dynamical friction time is assumed to depend on tlie cos- 
mic time and is shorter at higher redshift. This behaviour is ob- 
tained by using the Chandrasekhar formula with a multiplication 
factor of 5 X where t is the time in Gyr since the 
big-bang. This dependence is motivated by a more radial infall of 
satellite galaxies towards the centre of their group at high redshifts 
jHopkins et alj|201(]| : lweinmann et alj201 ih . 

• The SFR in merger induced bursts depends on the halo 
mass. In particular, it is l ower than the approximation used by 
lOe Lucia & Blaizo3 | |2007|) for haloes less massive than 3 x 
10^^ h~^MQ. The specific implementation is described below. 

Galaxies in our model are embedded within dark-matter sub- 
haloes as extracted from the A'^-body simulation. We assume that 
the galaxies merge following the merging-time of their host sub- 
haloes, with an additional time-delay. This delay is esti mated using 
the C handrasekhar formula for dynamical friction (see IWang et aU 
l2012h . Mergers trigger SF bursts, with an efficiency that depends 
on the mass ratio of the two galaxies. The total stellar mass formed 
in the merger induced burst is: 



0.56 



Am, 



T^T-cold 



0.56 ^ 

, mi 



T^T-cold 



Mo 



if Ah > Mo 



otherwise 



(1) 



where mi , m2 are the baryonic masses of the central and satel- 
lite galaxy respectively (including both stellai^ and cold gas mass), 
rricoid is the sum of the cold gas masses of the two galaxies, 
A'lh is the mass of the descendant subhalo, and Mo = 3 x 
10^^ /i^^Mq. For high mass galaxies, this recipe follows the re- 
sult s of hydrodynamic al simulations bv lMihos & Hemguisll ( 1 1994) 
and Cox et al. and was adopted by various SAMs (e.g . 

Somerville et al.ll200ll: ICroton et al.ll2006l iKhochfar & Silkll2009l: 



Neistein & Weinmann|[2oTol ; IKhochfar et al] 1201 ih . For low mass 



galaxie s this recipe is mo tivated by recent observations of SF effi- 
ciency jWang et al.l2012h . Note that mergers provide the only trig- 
ger for SF bursts in our model. The other mode of SF, by secular 
processes, describes the relatively slow conversion of cold gas into 
stars within disk galaxies. 

The specific model used here is able to fit the SF histo- 
ries of galaxies to a h igher accuracy than most other SAMs (see 
iHenrigues et al.ll2012l for a recent suc cess in this respect). As was 
shown in Neistein & WeinmanrjfcoiOl) . our mode s reproduce both 
the stellar mass function of galaxies up to z = 4, and the distribu- 
tion of SFRs at z = 0. In , Wang et al . (2 013) we further improved 
the model to reproduce the auto-correlation function of galaxies at 
z = 0. These results are unique in comparison to other SAMs, al- 
lowing us to use reliable SF histories of galaxies as a basis for the 
work presented here which focuses on BH evolution. The model 
fits observations to a level of 20-40 per cent, and is probably not 
the only possible galaxy formation model. In particular, there is a 
degeneracy between the amount of stars formed in mergers induced 
bursts, versus the amount of star formation due to secular processes. 

The duration of SF bursts has a negligible effect on the statis- 
tical properties of the model galaxies. This is because these events 
are rare, and hardly contribute to the total SFR density of the Uni- 
verse. However, in our current study, bursts are the only channel 
of BH growth, and their duration plays an important role in the 
model. We assume that the shape of SF bursts follows a Gaussian 
as a function of time, with the following parameterization: 



drua. 



AtTIs 



dt 



■ exp 



{t^tof 



2ai 



(2) 



Here cr(, determines the burst duration, t is the time in Gyr since the 
big-bang, and to is the time of the peak in the SF burst. We assume 
that to occurs 2ai, after the time the galaxies merge, to allow a 
smoothly rising peak of SFR (the burst of SF starts at the time of 
merging). We allow the burst to continue forming stars up to 4(Ti, 
after to, so the total duration of the burst is 6a i,. 

At higher redshifts galaxies have smaller radii, and time-scales 
are shorter (the halo dynamical time is proportional to the cosmic 
time, t). We therefore assume that at depends on t in the following 
way: 

t 



(3) 



where ao is a free constant. 

An additional important ingredient of the model is the defini- 
tion of the bulge mass. We assume that galaxies can grow a bulge 
according to two different channels. First, SF bursts that are trig- 
gered by mergers of any mass ratio are assumed to contribute their 
stellar content to the bulge (i.e., Ams.bmst from Eq. [T]is added 
to the bulge of the remnant galaxy). Second, once the mass ratio 
7712 /mi is larger than 0.3, the total stellar mass of both galaxies is 
moved to the bulge of the remnant galaxy. Therefore, the mass of 
stars within the bulge might be larger than the total amount of stars 
formed within bursts. 



2.2 Black holes 

Our SAM assumes that BHs grow only during merger induced star 
bursts. The masses of seed BHs in this model are extremely small 
in comparison to the mass added in the first merger event and make 
no difference to the accumulated mass of BHs at relatively low red- 
shifts. We have tested that using BH seeds of mass 10^ Mq does 
not change our model results significantly. 

When galaxies merge, we merge their corresponding BHs at 
the same time. Following each merger event, a burst of SF occurs 
according to Eas. lll3l We allow the remnant BH to grow in mass 
according to: 



dmsH 
dt 



-Oacc with probability Op 



otherwise 



(4) 



Here 77 is the fraction of mass that is transformed into radiation 
and aacc is a free parameter, corresponding to the efficiency of BH 
accretion with respect to the SF burst. 

In a ACDM universe, the fraction of galaxies with an active SF 
burst at any given time is high, especially due to the non-negligible 
burst duration at. However, obser vations show that th e number of 
AGNs is significantly lower (e.g. ICroom et al. I l2009h . Our model 
assumes that only a fraction ap of the merger events induce ac- 
cretion into the BH. In practice, at each time-step of the SAM 
(with a typical duration of 10 Myr) and for each BH we generate 
a random number, distributed uniformly between zero and unity. 
We then allow accretion only if this random number is smaller 
than Op. In practice, each merger event adds on average a mass 
of aaccOpAms, burst to the BH. 

The bolometric luminosity of the AGN is defined as: 



ri AmBH 2 
-t^AGN = z — — c 



(5) 



l-ri At 

where At is the length of the time-step within the SAM (~ 10 
Myr), AmBH is the total mass that is added to the BH within a 
time-step (i.e., the integration of Eq.|4j, and c is the speed of light. 
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The Eddington luminosity describes the limit at which radia- 
tion pressure balances the gravitational force of the BH 



1.5 



VlO^Me 



10 ergs 



(6) 



where the factor 1.5 on the left is derived for solar metallicity gas. 
In our model, the accretion into the BH is determined by the SF 
burst, and might be above the Eddington limit. In addition, it is rea- 
sonable to assume that this limit will vary between different BHs. 
In order to use a practical accretion limit, we define Asdd to be 
a log-normal random variable (i.e. its log value is normally dis- 
tributed) with a mean of AecIcI.o and a standard deviation of (TEdd 
(in log). For each accretion event, we randomly choose a value for 
AEdd and do not allow Lagn to exceed LEddAEdd. We will refer 
to LsddAEdd as the Eddington threshold, since it is being used to 
actively limit the accretion onto BHs. We note that we do not use 
the Eddington threshold for the first accretion event (i.e. when the 
BH mass equals zero). 

At low accretion rates we use the properties of 'advection 
dominated accretion flows' (ADAFs). As in previous studies (e.g. 
iFanidakis et alj2012l) we model the ADAF limit by 

I/AGN 



»7ADAF — 



V 



(7) 



ifidd QADAF 

which is valid only for Lagn / Lsdd < Qadaf. By definition, 
?7adaf equals rj when I/AGN/iEdd = oadaf- In this work we 
fix oadaf to a value of 0.01. 

Some of the results presented here make use of the luminosity 
of AGNs in the B band. For the conversion between Lagn and B 
band we use the bolometric correction fromlMarconi et aL](|2004t): 

log {i^bL^^ / Lagn) = -0.8+0.067Li-0.017i?+0.0023L? .(8) 

where Li = log(LAGN/Lo) — 12. The absolute magnitude in the 
AB system is then given by 

A/s = -11.32 -2.5 log (!^sL,.b/10'"' ergs"') , (9) 

where B refers to a rest wavelength of 4400A. 

To summariz e, our SAM includes all the ingredients from 
IWang et al.l ( [2012h related to galaxies, with the following additional 
parameters related to the evolution of BHs: 77, aacc, ctp, Asdd.o, 
and (JEdd. These parameters are presented in Table [Tjand their best 
chosen values are justified in sectio n [3] The ingredients used here 
are similar to previous SAMs (e.g. jMalbon, et al ,[ |2Q0 ?b , although 
we do not include grow th modes that are not due to bursts (e.g. 
iHirschmann et al 

Finally, we note that the model is limited both in terms of the 
minimum BH mass that is properly resolved, and in terms of the 
limited volume of our simulation box. These limitations originate 
from the limited dynamical range of the Millennium simulation that 
is being used here, and will be discussed below. 



3 RESULTS 

3.1 Observational constraints and best model parameters 

In this work we treat the galaxies as priors and do not change the 
parameters that affect their evolution, except for the value of ctq, 
which has a negligible effect on the properties of galaxies. Conse- 
quently, there are only six free parameters in the model, as listed in 
Table[T] However, the low density of AGNs in the Universe forces 
us to test each set of model parameters by using the full Millen- 
nium simulation, spending a few CPU hours on each run. In order 
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Figure 1. The mass functio n of BHs. The shaded region represents the ob- 
servational prediction from Ishankar et al.| (|2009j) based on the congelation 
between bulge and BH masses at low redshift. Different line types corre- 
spond to the model mass function at various redshifts as indicated. The grey 
dashed line shows the minimum BH mass that is reliably reproduced by the 
model (due to a minimum subhalo mass of 1.72 X 10^° h-iMo). 



to achieve fast tuning, while still running our model on the full 
simulation box, we first save all the model results that are related 
to the evolution of galaxies. We then run the BH ingredients only, 
thus saving a large amount of computational time. Our tuning pro- 
cedure can evolve BHs over the full Millennium simulation in only 
three minutes (using one processor), allowing us to systematically 
explore a large region of the parameter space. 

We use three different types of observations to constrain our 
model parameters: the luminosity functions of AGNs at < z < 4, 
the mass function of BHs at z = 0, and the relation between SFRs 
and AGN luminosities. While we further compare the model to 
other observations in the following sections, these additional ob- 
servations were not used for tuning the model. 



3.1.1 The mass function of BHs and the luminosity functions of 
AGNs 

We start by tuning the model parameters against two common ob- 
servations: the mass function of BHs at redshift zero, and the lu- 
minosity function of AGNs at z < 4. We apply a tuning proce- 
dure that systematically scans a range of the parameters aacc, ctp, 
AEdd.Oi cTEdd according to the range listed in Table[T] The value of 
(70 can significantly affect the luminosity function of AGNs. How- 
ever, it is fixed here due to its effect on the relation between SFRs 
and AGN luminosities, that will be discussed in the next section. 
Using the prior value of cro = 0.13 Gyr, we could not find a model 
that matches the data to a satisfying accuracy, while using constant 
values for all the parameters. We note that larger values of ctq allow 
for better models. Consequently, we choose to add a dependence 
on time for aacc- It should be noted that a similar result as shown 
here could arise from varying the parameter ap with time. We will 
discuss below the degeneracy in choosing the best solution, and its 
relation with the SFR values of our model galaxies. Our best model 
uses rj = 0.05, consistent with a value of 0.057 obtained for non- 
rotating BHs. 

In Fig. [T] we show that the mass function of BHs from the 
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Table 1. The free parameters used to model BHs and AGNs. 



Parameter range best value description 





0.1-0.5 Gyr 


0.13 Gyr 


Burst duration, Eq. |3] 


V 


0.05-0.3 


0.05 


Radiative efficiency, Eq.|4] 


^acc 


0.001-0.3 


0.185(tf - 1.68*1 - 


- 0.75) Accretion efficiency, Eq.|4] ti = */13.6 Gyr 


Ctp 


0.001-0.1 


0.04 


Fraction of mergers that result in AGN activity, Eq.|4] 


-^Edd.O 


0.1-3 


1.0 


Eddington threshold (mean) 


0"Edd 


0-1 


0.7 


Eddington threshold (scatter) 



X -9 

CP 

CO —5 

I 

a -7 



"-9 

bJO -5 

o 

^ -7 
-9 



z=0.4-0.6r* 


z=0.68-1 .06^* 


z=1 .06-1.44 *^ , 








z=1 .44-1 .82 
1 1 1^ 


z=1. 82-2.20 


z=2.20-2.60 % ^ 








z=3.25 ^ , 


s 
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z=4.25 ■ 
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Figure 2. The luminosity functions of AGNs. Each panel shows one redshift bin as indicated. Observational results from lCroom et aljj2009l) and lRichards et al] 
)200d) are plotted in diamonds and squares symbols respectively. Following Eq. 2 and 3 from lCroom et al.l )2009h . we add 0.587 and 0.817 magnitudes to Mg 
and Mi respectively, in order to convert them to Mg (note that in these papers Mg and Mi are defined for rest wavelengths corresponding to z = 2 and hence 
require correction factors to compare with our computed B band at a wavelength of 4400A). Model results are plotted in solid lines. 



model at z = agrees with the predi c tion f rom IShankar et al.l 
( I2OO9I) . The results from IShankar et all (|2009|) are based on the 
mass function of galaxies, combined with the correlation between 
the bulge and the BH mass. It is therefore not a direct observa- 
tion, and should not be considered as a strong constraint. However, 
the mass function points to important features of the model and we 
choose to use it as a constraint even thoug h the actual values are no t 
certain. For example, the recent work by iMcConnell & Mai ( I2OI2I) 
indicates that the ratio between the mass of bulges and BHs could 
reach a level of ~ 100 at the high-mass end. Fig.[T]also shows the 
model mass functions at higher redshifts, indicating that the num- 
ber of low mass BHs grows fast before z ~ 1 and changes only 
slightly at z < 1. 



As seen in Fig. [Tj the computed mass function declines at 
logmBH/M© < 7.7. We have tested a lower resolution model 
by allowing BHs to grow only within subhaloes more massive 
than 10" h^^MQ. This mimics a low-resolution simulation with 
respect to the actual minimum mass of the simulation (1.72 x 
10^" /i~^Mq). Within this resolution test, the bend in the mass 
function occurs at a larger mass, indicating that our resolution limit 
is indeed the reason for the bend at log rriBH /Mq < 7.7. Our lim- 



ited resolution has a negligible effect on other results shown in this 
work. 



Models based on high resolution me rger trees could possi - 
bly solve the above resolution problem (e.g. lFanidakis et al]|2012h . 
Since it is difficult to increase the mass resolution within cos- 
mological simulations, this issue could be addressed by using 
Monte-Carlo algorithms for generating merger-trees of haloes (e.g. 



Somerville & Kolatlll999l:|Parkinson et alj2008l:lNeistein & Dekeil 



20081) . However, the merger-trees used here are based on subhaloes 
and are more accurate than the theoretical merger trees based on 
haloes. We plan to check the use of Monte-Carlo merger trees in a 
future study. 

The luminosity functions (LFs) of AGN s are shown in Fig. 2| 
where points are obs ervations adopted from ICroom et al.l (12009 ). 
iRichards et all ( |2006|) and solid lines are the model results. The 
theoretical LFs computed by the SAM are in broad agreement with 
the observations. In case we use constant otacc the LFs at z < 1 are 
higher than the observed ones. 

It should be noted that our simulation volume is smaller than 
the observed volume at z > 2, giving rise to some of the deviations 
seen at high luminosities at this redshift range. Furthermore, we do 
not take into account various effects that could change the model at 
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Figure 3. The average SF R for host galaxies o f a given AGN luminosity. 
Observational results from iRosario et an i l2012l) are shown in grey dashed 
lines, where thinner lines con'espond to higher redshifts. The observed en'or 
bars are typically at the level of ±0.2 and ±0.5 dex, for SFR and Lagn 
respectively. Model results using the same averaging scheme and the same 
redshift bins are plotted in solid lines with redshifts as indicated. For refer- 
ence, we plot in dotted line the relation SFR oc i^GN (''^^ P^"*- refen'ed 
to in the paper as the 'rates diagram'). 



the level of ~30 per cent. For example, it is well known that some 
AGNs are obscured and would not be accounted for in the observed 
LFs. We have tested that using the slim-disk approximation for su- 
per Eddington accretion rates (e.g. Fanidakis et al. 2012) does not 
change the model results. Lastly, re-tuning the ingredients related 
to galaxy-formation might further improve the model results. We 
neglect all these effects in order to keep our model as simple as 
possible. 

Most of the free parameters in our model have a simple effect 
on the LF. Since the SF bursts of galaxies are not modified, Oacc 
behaves as a constant multiplication factor for the luminosity of all 
AGNs at a given redshift. Changing aacc thus induces a lateral shift 
in the luminosity function. The value of ap changes the number of 
bursts leading to BH activity, and therefore shifts the luminosity 
function along the Y-axis. Both aacc and Op are degenerated at 
some level, yielding similar results for different combinations with 
the same value of aaccCtp- The value of rj is similar in its effect to 
Uacc- In addition, it changes the relations between mass growth (i.e. 
the mass function of BHs) and luminosity. The parameters of the 
Eddington threshold, AEdd.o and a-^dd, change mostly the high-end 
part of the luminosity functions, and also the mass function of BHs. 
These two parameters do not have a simple effect on the model 
results. When neglecting the Eddington threshold, the luminosity 
of AGNs only depends on the SFR within bursts, and not on the 
BH mass. In this simplified case, the LF of AGNs has the same 
shape as the number-density of SFRs in the model. 



3.1.2 Star-formation rate and AGN luminosity 

In Fig.[3]we show the observed results from lRosario et al.](l2012h in 
grey dashed lines. Values along the X-axis correspond to the bolo- 
metric luminosity of X-ray selected AGNs. The original Y-axis of 
this plot corresponds to the rest-frame luminosity at a wavelength 
of 60 fim, stacking all galaxies that host AGNs with similar Lagn. 
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Figure 4. The rates diagram: values of SFR versus Lagn at ; = 0. Cir- 
cles represent the SFRs and AGN luminosities of the model galaxies. Cross 
and square symbols show the level of SFR for the same objects, when con- 
sidering contributions from SF bursts only (symbols are plotted for 20 per 
cent of the objects). We use the variable x = [t — to)/ai,, where x < —1 
corresponds to objects in the early stages of the burst. The thick solid line 
shows the average SFR value for a given Lagn. was computed in Fig.|3] 
(i.e., based on the total SFRs). The dashed line corresponds to the average 
value of Lagn . for ^ given value of total SFR. 

In order to transform the IR luminosity into SFR, we use the rela- 
tion L/h = 3.8x 10'''^SFR (erg s^^) and Lbom™ — 0.5I//ii, where 
SFRs are given in units of Mq yr^^. The corresponding SFR val- 
ues are thus population averages for various redshift groups and do 
not represent individual objects. 

Interestingly, as discussed in iRosario et all | |2012|) , for z ^ 1 
there is hardly a correlation between the AGN luminosities and 
SFRs. At z < 1 the correlation between SFRs and AGN luminosi- 
ties only exist for Lagn > iC'' erg s^^. These findings seem 
surprising, since we expect that AGN growth modes will be cou- 
pled to the properties of their host galaxies. A lack of correlation 
might indicate that most AGNs are triggered by processes that are 
not related to the dominant mode of SF. In addition, it does not 
seem obvious that the congelation for bright AGNs only exist at low 
redshifts. In what follows, we term this plot as the 'rates diagram', 
since it shows the relation between BH accretion rates, and SF rates 
of their host galaxies. 

An explanation for the rates diagram should take into account 
the various time-scales that are buried within this relation. First, 
the observations of SFRs are based on indicators that last ~150 
Myr (this is a combination of both UV light emission from young 
stars, and dust heating that yields IR luminosity). Hence the ob- 
served SFR values reflect some of the recent history of the galaxy. 
To follow the same effect in the model, we compute the average 
SFR for each galaxy over time, taking into account the last 150 
Myr before the output snapshot. Second, AGN luminosities repre- 
sent an instantaneous emission, assumed by the model to last for 
~ 10 Myr. Third, the duration of both the burst of SF and AGN 
activity is 60-6 « 0.06t (see Eqs.|2]&[3]l. 

As a result of the SF time-scale, the observed value of the IR 
luminosity is affected by the secular SF that occurred before the 
burst. This luminosity is uncoiTelated with the merger event and 
with the AGN luminosity. High luminosity AGNs are those that 
are observed at a time that is close to the peak of their accretion 
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Figure 5. The temporal location of AGNs with respect to the peak lumi- 
nosity of the burst. For each value of SFR and Lagn we plot the median 
value oi X = {t - to) /at (see Eq.|2). All AGNs are selected at redshift 
zero. The dashed line follows the relation SFR oc i^GN- AGNs that are 
starting their accretion episode lie on the right part of the diagram. 



event (to from Eq.|2]l- In these cases the SF burst is significant, and 
consists of most of the observed SFR value. This explains why the 
congelation between SFRs and AGN luminosities is only seen at 
high AGN luminosities. At z = 2, the secular SFR for star-forming 
galaxies is high, reaching the same level of SFR due to bursts in the 
most luminous AGNs. This is the reason for the lack of correlation 
between SFR and AGN luminosity at these redshifts. 

At a given observed epoch t, different bursts have different 
values of to and Ams^burst only, as all other parameters do not vary 
between different objects. Since the time-scale for averaging the 
AGN luminosity is short, objects with the same value of 

^^s, burst 

and \t — to\ will have the same values of Lagn- However, these 
objects will have very different SFRs because those with t > to 
will include the peak of SF within their SFRs, while for t < to the 
SFR will mostly include the secular SF. This effect is more severe 
for larger values of \t — to\, and once the burst duration is smaller 
than the time-scale of SFR, contributing significantly to the scatter 
at low Lagn- 

In Fig.|4]we demonstrate these effects on our model galaxies 
at z = 0. We first plot individual SFRs and AGN luminosities from 
the model, and show how the average SFR is flat at Lagn < 10''* 
erg s~^. We then plot for the same objects the SFRs that are cal- 
culated based on bursts only, taking out the contribution from the 
secular SF mode. The SFR due to the bursts is more correlated with 
the AGN luminosity. The remaining scatter between the bursty SFR 
and Lagn originates from the fact that AGNs can be selected both 
before and after the peak (i.e. to can occur before or after the ob- 
servation). This changes the contribution to SFR from the burst it- 
self, although the AGN luminosity is left unchanged. In Fig.|4]we 
show in different symbols AGNs that are selected before the peak 
of the burst, with (t — to)/o"6 < —1- For these objects, the correla- 
tion between SFR and AGN luminosity is strong, as expected. The 
strength of this last effect is smaller when at is higher (even when 
using (To — 0.2 instead of the current value of 0.13, the scatter due 
to different values of x becomes less dominant). 

Our model suggests that only ctq and aacc affect the rates di- 
agram. Changing the value of ctq results in modifying the average 



level of SFR for low Lagn- On the other hand, Uacc shifts the 
diagram along the X-axis. Unlike the behaviour of the luminosity 
function, there is no simple degeneracy between these parameters 
here. The value of cro = 0.13 Gyr that is chosen for the best model 
is obtained by matching the observed rates diagram. 

We have tested how other standard SF modes affect the rates 
diagram and found that the Lagn-SFR correlation is very sensi- 
tive to secular modes. For example, when using accretion into BHs 
that is proportional to the secular SFR, there are no regions with 
uncorrelated SFR and AGN luminosity. This is apparent even in 
cases where the BH growth in the secular mode is only 10" '' of the 
assumed SFR. As indicated by Fig.|4l the rates diagram originates 
from two effects. First, there is hardly a correlation between the SF 
within the burst and the secular SF. Second, the episodic nature of 
the burst together with the integration time of the SFR indicator. 
These two effects do not exist in a secular mode that is based on 
continuous BH accretion mode. However, different types of secu- 
lar accretion, e.g. episodic events that were not tested by us, might 
show better agreement with the observed rates diagram. 

The rates diagram as shown here is defined as the average SFR 
at a given AGN luminosity, and is not affected by the number of 
AGNs of different luminosities. However, once we compute the av- 
erage AGN luminosity as a function of SFR (dashed line in Fig.|4]l, 
we see a different trend, in which SFR and Lagn are always cor- 
related. This way of averaging hides the flat part of the diagram, 
occupied by low luminosity AGNs, and is more dependent on the 
number of low-luminosity AGNs (i.e. on sample selection). Since 
the flat part of the rates diagram is highly restricting the models, 
the average Lagn is less restrictive when comparing models to ob- 
servations. 

In Fig.[5]we plot for each value of SFR and AGN luminosity 
the median value of x = {t — to) /ai,. According to Eq.|2] x is the 
location within the Gaussian of the burst, where negative values 
(down to -2) correspond to the burst starting point, and positive 
values (up to 4) con'espond to the late episode of the burst. Each 
burst starts on the left part of the rates diagram, with small AGN 
luminosity that is due to a value of a; = —2. As a; increases towards 
0, both the SFR and Lagn rise. Later on, when a; > Lagn goes 
down to its initial level, while the SFR stays a bit higher, due to the 
integration over 150 Myr. As a result, the right part of the diagram 
is mostly populated with a; < — 1, and the higher values of Lagn 
are related to a: ~ 0. This effect modifies the correlation of SFR oc 
Lagn (that is used by the m odel, Eq.|4]l such that SFR oc Lagn> 
as observed bv lNetzerl l l2009l) . 

3.2 Merger fraction 

The main argument against a model of BH growth that is based 
only on mergers is the morphology of AGN hosts. As discussed in 
the Introduction, numerous works show that most AGNs reside in 
normal disk galaxies, with no sign of interaction. Here we test the 
observational constraint in this respect more closely. 

In Fig.[6]we show the fraction of galaxies that host AGNs and 
are participating in major merger events, out of all galaxies that host 
AGNs. Major merger events are defined as those with mass ratio 
that is bigger than 0.3. Note that in the SAM, if a galaxy have more 
than one merger event, the mass ratio is determined by the event 
with the largest contribution to the SFR. A galaxy is considered to 
be a part of a merger event throughout the burst. 

Fig.[6|shows tha t the results of the model are roughly consis- 
tent with those oflTreister et al. (2012.) . although they deviate from 
the analysis of lKocevski et alj j2012h . The model predicts that ma- 
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Figure 6. The fraction of galaxies that host AGNs and are going through 
major mergers (mass ratio is bigger than 0.3) out of all galaxies that host 
AGNs. Thick lines show the model results at various redshifts as indicated, 
after selecting all galaxies that host BHs with masses log mBH/M(7) > 
7.7. T he thin dashed line shows the fitting function from iTreister et al.l 
obtained from a compilation of various studies at < z < 1. 
The grey boxy line shows the results from lKocevski et al] )2012h . derived at 
z ~ 2. 
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Figure 7. The fraction of galaxies undergoing a merger event with a mass 
ratio r. dP/dr is the probability distiibution of galaxies of different r values. 
All galaxies are selected at z = 2.1 with stellar mass as indicated, in order 
to avoid resolution issues. The solid line corresponds to all the galaxies 
within the model. The dashed line shows the fraction out of all the galaxies 
th at host an AGN with 1 0^^ < Lagn < 10'''' erg s~l, as were selected 
bv lKocevski et alj <2012h . 



jor mergers are associated with more luminous AGNs. There are a 
few limitations in this comparison that should be mentioned. First, 
we assume that the mass ratio for major mergers is 0.3, which is 
somewhat arbitrary. Second, we assume that the time-scale for de- 
tecting the morphological features within the galaxy is similar to 
the assumed duration of the burst. Third, we do not allow for a 
time-delay between the change in morphology and the onset of ac- 



cretion into the BH. All these issues complicate the interpretation 
of the results. 

A different way to examine this issue is to compare galaxies 
that host AGNs to a control sample of galax ies with the same stella r 
mass, that do not host AGNs. According to lKocevski et alj ( l2012h . 
these two populations show similar fractions of galaxies that are 
detected as going through merger events. In Fig. |7]we make the 
same comparison with our model galaxies. This plot shows the 
fraction of galaxies going through a merger event of mass ratio 
r, after selecting only galaxies at z = 2, and with stellar mass 
10.5 < log(mi,/MQ) < 11. The AGN sample includes galax- 
ies hosting an A GN with 10"'^ < Lagn < iC"' erg s"^, similarly 
to the sample o f lKocevski et"ai] j2012h . This result shows that the 
fraction of galaxies going through a merger event is larger by a 
factor of ~ 2 for galaxies that host AGNs. Unlike the comparison 
made in Fig. (6] here the uncertainty is due to the existence of a 
time-delay between changes in morphology and the BH accretion 
phase. 



3.3 Bulge and black hole mass 

In Fig. [8] we show the relation between bulge and BH mass in our 
model. The ratio of the two masses depends on various factors. 
First, bulges can grow in major merger events by using the stellar 
mass from the disks of both progenitor galaxies. This effect will 
allow bulges to gain more mass than the total SFR during mergers, 
more than what is used for feeding the BH. In addition, the scatter 
between the bulge and BH masses should be a result of ap, which 
allows for only a part of the burst events to feed the central BH. 

Our model agrees wit h the results of both iHaring & Rixl 

( l2004b and lSani et al. 1 1(201 ij) atredshift zero. However, the recent 
study by iMcConnell & M j {2Ql% shows that at a bulge mass of 
lO^'^ M0 the BH mass is higher than what we get here, reaching 



10 Mq. Our model BHs are less mas sive at the high mass end 
because we use the mass function from IShankar et al.l ( |2009|) for 
calibrating the mass of BHs in our model. In addition, the model 
predicts a smaller scatter in the mass relation for more massive ob- 
jects in comparison to observations. 



3.4 The distribution of BH masses for a given AGN 
luminosity 

The distribution of BH masses, for a given value of Lagn are 
shown in Fig. |9] The model predicts that most AGNs accretion is 
below the Eddington limit, especially at low redshift. At higher red- 
shifts and high Lagn the AGNs accretion is closer to the Edding- 
ton limit. The scatter we use for the Eddington threshold allow for 
some objects to slightly exceed the formal Eddington limit, mostly 
at high redshift. 



3.5 Specific SFR in AGN hosts 

Interestingly, our merger scenario naturally predicts that AGNs 
are hosted by massive star-forming galaxies. This can be seen in 
Fig. Tol Observations of this type were done bv lSalimetai1 ( l2007l) 
and Bongiomo et al.l ( I2OI2I) using different SFR indicators. This 
diagram is an extension of the rate diagram shown in Fig. |3] as it 
shows the full distribution of specific SFR, for each value of stellar 
mass. Results at higher redshifts are similar to what is shown here 
for z = 0. 
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Figure 8. The relation between the masses of bulge s and BHs at z = . 
Sy mbols correspond to the observational results from iHaring & Rixl | |2004) 
and lSani et alj i201 lb . The thick solid line shows the median BH mass for a 
given bulge mass, taken from the model galaxies. Dashed lines show the 5 
and 95 per cent levels of the BH distribution for each bulge mass. 



Figure 10. The distribution of specific SFR (SSFR) versus stellar mass 
at z = 0. Thick contours represent the distribution of all the model 
galaxies, while thin lines show the distribution of galaxies that host 
AGNs with Lagn > 10*^ erg s^^. Contour levels ai'e at logP = 



-2.9, 



-1.4, —0.9, —0.4, 0.1 where the integral of P equals unity. 



4 SUMMARY AND CONCLUSIONS 
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Figure 9. The distribution of BH masses as a function of AGN luminosities 
at various redshifts. The thick solid lines are the median values of the model 
The dashed lines show the 67 and 95 per cent levels of the distribution. The 
thick straight line in each panel shows the Eddington luminosity, Eq.|6] 



3.6 BH mass and accretion rate in the SFR-Lagn plane 

In Fig.[TT|we show how other properties of AGNs depend on both 
SFRs and AGN luminosities at z = 0. We show that the model 
predicts a higher values of I/AGN/iEdd at higher Lagn, with 
small dependence on SFRs. This is in spite of the lower values of 
X = {t — to) /at seen in Fig. [5] at high Lagn and low SFRs. As 
a result, the Eddington ratio does not depend on x, the time within 
the burst. We have further tested that the Eddington threshold does 
not play a role at shaping the rates diagram at z = 0. 

In the lower panel of Fig. [TT]we show the median BH mass 
at each value of SFR and Lagn. The scatter in masses around 
Lagn < 10^'* erg s^^ depends on the SFR but other regions of 
the diagram include more uniform distribution of BH masses. 



This work presents a semi-analytic model (SAM) for the formation 
of galaxies and black holes (BHs) within a ACDM Universe. Our 
model is built on a galaxy-formation model that matches various 
observations of galaxies. These include the stellar mass function 
at < z < 4, the distribution of star-formation rates (SFRs) at 
z = 0, and the auto-correlation function of galaxies at z = . The 
results were developed in a previous study (Wang et al.ll2012l) . and 
are adopted here with no further modifications to the properties of 
galaxies. When adding BHs to the model, we are interested in the 
properties of BHs and active galactic nuclei (AGNs), and how they 
coiTelate with the properties of their host galaxies. We do no take 
into account the possible feedback from AGNs on the physics of 
gas within galaxies. 

Our model assumes that BHs grow only during star-formation 
bursts, with a probability Op that is the same for all galaxies at all 
redshifts. The bursts in our model are only due to merger events; 
we do not model bursts from secular processes such as disk in- 
stability. When an event of BH accretion occurs, we assume that 
the gas mass that reaches the BH equals a fraction Oacc of the 
stars made in the burst at the same time. We use two additional 
parameters to describe the threshold on accretion due to the Ed- 
dington luminosity, in order to limit the accretion in case the burst 
is large and the BH mass is small. This model for the formation 
and evolution of BHs is therefore very simple, also in compari- 
son to previous S AMs ( Kauffmann & Haehnelt 2000; Croton e t al] 

2OO7F 



Bo wer et al]|200 
Somerville et al. | I2OO8I ; iBonoh et al 



2006; 



Malbon et al. 2007; Monaco et al. 



2009,; .Fanidakis et al. 



2011 



Fanidakis et al.ll2012l ; Hirschmann et al.ll2012ir 

The observational constraints used here are the mass function 
of BHs at redshift zero dShankar et al.il2009l) the luminosity func- 
tion of AGNs at z < 4 jRichards etal.1120061 ; ICroom et alboogh. 



and t h e relation between S FRs and AGN luminosities ( Shao et al] 
I2OIOI ; iRosario et al.ll2oTi) . We specifically study the latter con- 
straint, as it was n ot tested by previous models. According to 
iRosario et alj ( |2012|) . AGNs at low redshift are related to a con- 
stant level of SFR within their host galaxies for Lagn < lO*** erg 
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Figure 11. SFR versus Lagn at z = 0. The upper panel shows the median 
Eddington ratio, log Lagn / LBdd. at each bin. The lower panel represents 
the median BH mass in units of log Mq. Here we set BH masses to a min- 
imum value of 10^ Mq . 



event, and the amount of star-formation within the burst, which de- 
pends on the properties of both merging galaxies. We thus argue 
that the rates diagram favours a mode of BH accretion that is due 
to mergers. 

The rates diagram is used here for constraining the duration 
of the bursts in our model and for constraining the BH accretion 
factor, ctacc - The diagram is very sensitive to secular modes of BH 
growth (i.e. modes that do not depend on merger events). Even a 
small level of BH accretion from the secular mode, which adds a 
negligible mass to BHs, can modify the rates diagram considerably. 
Although we do not model such scenarios here, more constraints 
could be obt ained by measuring the SFR in the central regions of 
galaxies (e.g. lStem & Laoiil2012h . 

Visual classification of AGN hosts indicates that most AGNs 
are located within normal disk galaxies, w i th no signs of merg- 
ers jGabor et al.ll2009t ICistemas et al1|201 ll: iKocevsld et al.ll2012l : 



ISchawinski et al 



201^ 



Treister et al.l 120121) . We show in Figs. |6] 



&|7]that a similar result is obtained here, using a model with BH 
growth coming only from mergers. This is because most BH ac- 
cretion is associated with minor merger events, which are difficult 
to detect observationally. However, our model shows that a typical 
galaxy has roughly half the chance to experience a merger event in 
comparison to a galaxy of the same stellar mass that hosts an AGN. 
We do not test in this work the possibility of having a time-delay 
between the merger event and the episode of BH accretion. Such a 
time-delay will further decrease the amount of AGN hosts that are 
showing morphological signs of interactions. 

The model parameters used here were tuned to fit the ob- 
served mass function of BHs at z = 0, and the luminosity func- 
tion of AGNs at z < 4. In addition, our final model agrees with 
the observed correlation between the mass of bulges and BHs 
( iHaring & Rixl2004ISani et al.l201 ihlMcConnell & Mj2012h . We 
show in Fig. [To] that our model predicts that AGN hosts are rela- 
tively massive, and star-forming galaxies. 

Finally, we provide predictions for the mass of BHs, and their 
Eddington ratios, within the rates diagram, as a function of both 
SFR and Lagn ■ These predictions could be used to test the model 
against new observations. 



s~^. At higher AGN luminosity, where Lagn is larger than the lu- 
minosity of S FR, a power-l aw type correlation, SFR oc Lagn is 
seen (see also lNetze3l2009l) . At z ~ 2, AGNs of all luminosities 
are hosted within galaxies of the same average SFR, at a level of 
~ 30 M0 yr~^. We term these results as the 'rates diagram', since 
it shows the relation between the star-formation rate, and the rate 
of BH growth. 

As was demonstrated in Figs.[3]&|4] our model fits well the 
observed rates diagram. A necessary condition for the agreement 
is that we adopt the correct time scales for both SF and AGN ac- 
tivity. We the r efore assume that the luminosity at 60/im used by 
iRosario et al. I l l2012l) represent an averaging of 150 Myr over the 
SFR of each galaxy, and that the AGN luminosity is instantaneous 
(i.e. using the minimum time-step of our model, ~ 10 Myr). When 
using these time scales, we naturally reproduce the observed rates 
diagram. In our model, the lack of correlation between SFR and 
AGN luminosities is due to two main reasons: First, secular SF in 
galaxies are not associated with BH activity. Second, BH accretion 
before and after the peak in SFR could have the same AGN lumi- 
nosity, but different SFR. A key assumption in the model is the lack 
of correlation between the properties of a galaxy prior to the merger 
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